Braneworld tensor anisotropies in the CMB 
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Cosmic microwave background (CMB) observations provide in principle a high-precision test of 
models which are motivated by M theory. We set out the framework of a program to compute 
the tensor anisotropies in the CMB that are generated in braneworld models. In the simplest 
approximation, we show the braneworld imprint as a correction to the power spectra for standard 
temperature and polarization anisotropies. 
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The early universe provides a testing ground for theories of gravity. The standard cosmological model, based on 
general relativity with an inflationary era, is very effective in accounting for a broad range of observed features of the 
universe. However, the lack of a consistent theoretical framework for inflation, together with the ongoing puzzles on 
the nature of dark matter and dark energy, indicate that cosmology may be probing the limits of validity of general 
relativity. 

M theory is considered to be a promising potential path to quantum gravity. As such, it is an important candi- 
date for cosmological testing. In the absence of a sufficiently general M-theoretic model of cosmology, we can use 
f^i , phenomenological models that share some of the key features of M theory, including branes. In brane cosmology, the 
observable universe is a 1+3-dimensional "brane" surface moving in a higher-dimensional "bulk" spacetime. Standard- 
model fields are confined to the brane, while gravity propagates in the bulk. The simplest, and yet sufficiently general, 
phenomenological braneworld models are those based on the Randall-Sundrum II scenario [1]. These models have 
the additional advantage that they provide a framework for investigating aspects of holography and the AdS/CFT 
■ correspondence. 

In the generalized RSII models, analyzed via an elegant geometrical approach in Ref. [2], the bulk is a 1+4- 
dimensional spacetime, with non-compact extra spatial dimension. What prevents gravity from 'leaking' into the 
extra dimension at low energies is the negative bulk cosmological constant A5 = —6/£ 2 , where £ is a curvature scale 
of the bulk. In the weak-field static limit, null results in tests for deviations from Newton's law impose the limit 
mm. The negative A5 is offset by the positive brane tension A, which defines the energy scale dividing low from 
high energies. The limit I < 1 mm implies A > (100 GeV) 4 , and the effective cosmological constant on the brane is 
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A = i(A 5 + K 2 A), (1) 

where k 2 = 8ttG = 8tt/-MJ, and M4 ~ 10 19 GeV is the effective Planck scale on the brane. A further intriguing 
feature of the braneworld scenario is that, because of the large extra dimensions, the fundamental energy scale of 
gravity can be dramatically lower than the effective Planck scale on the brane - as low as ~ TeV in some scenarios. 
In generalized RSII models, the fundamental scale is higher, M5 > 10 5 TeV, and is related to M4 via Mf = M\jl. 

At energies well above the brane tension A, gravity becomes 5-dimensional and significant corrections to general 
relativity occur. There are also corrections that can operate at low energies, mediated by bulk graviton or Kaluza-Klein 
(KK) modes. Both types of correction play an important role in tensor perturbations. 

The background cosmological dynamics of a Friedmann brane in Schwarzschild-Anti de Sitter (AdS) bulk are well 
understood [3] , including the high-energy modifications to inflation [4] . High-energy inflation on the brane generates 
a zero-mode (4D graviton mode) of tensor perturbations, and stretches it to super- Hubble scales. This zero-mode has 
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the same qualitative features as in general relativity, remaining frozen at constant amplitude while beyond the Hubble 
horizon, but the overall amplitude is higher [5]. The massive KK modes (5D graviton modes) remain in the vacuum 
state during slow-roll inflation. The evolution of the super-Hubble zero mode is the same as in general relativity, so 
that high-energy braneworld effects in the early universe serve only to re-scale the amplitude. However, when the zero 
mode re-enters the Hubble horizon, massive KK modes can be excited. Qualitative arguments [6, 7] indicate that this 
is a very small effect, but it remains to be properly quantified, so that the signature on the CMB may be calculated, 
and constraints may be imposed on the braneworld parameters. 

We develop here a formalism to compute the tensor anisotropies in the CMB, which incorporates the early-universe 
high-energy braneworld effects, and we carefully delineate what is known on the brane from what is required from bulk 
equations. Once the 5D solutions are provided, our formalism, with its modified CMB code (based on CAMB [8, 9]), 
is able to compute these anisotropies. We illustrate this by using a simple approximation to the 5D effects (cf. the 
analysis of the braneworld scalar Sachs- Wolfe effect in Ref. [10]). 



II. BRANEWORLD DYNAMICS AND TENSOR PERTURBATIONS 



There has been an explosion of interest in the theory of cosmological perturbations in braneworlds and their 
implications for observational cosmology (see e.g. Refs. [11]). 

In this paper, we follow the same general formalism previously developed to analyze braneworld scalar anisotropies 
in the CMB [12]. This is based on the 1+3-covariant perturbation theory for CMB anisotropies [13], generalized to 
incorporate braneworld effects. (See Ref. [12] for further details and definitions.) 

The field equations induced on the brane are 

K 2 

Gab = -A 5afc + K 2 T ab + 6 — Sab - Sab , (2) 

where S a b ~ {Tab) 2 carries high-energy corrections, and is negligible for matter energy densities p« A, while £ a b, 
the projection of the bulk Weyl tensor on the brane, carries corrections from KK or 5D graviton effects. From the 
brane-observer viewpoint, the energy-momentum corrections in S a b are local, whereas the KK corrections in S a b are 
nonlocal [15, 16], since they incorporate 5D gravity wave modes. These nonlocal corrections to the Einstein equations 
on the brane cannot be determined purely from data on the brane, and so the induced field equations are not a closed 
system; one needs to supplement them by 5D equations governing £ a b- 

The trace free E a b contributes an effective energy density p* , pressure p*/3, momentum density q* and anisotropic 
stress 7r* b on the brane, which collectively incorporate the spin-0, spin-1 and spin-2 modes of the 5D graviton, and 
through which the brane "feels" the bulk gravitational field. Then the braneworld corrections can conveniently be 
consolidated into an effective total energy density, pressure, momentum density and anisotropic stress: 

P° ff = pfl + £ + ^ . (3) 



2A p 

if = g (i + £)+?: (5) 

off ( , P + 3P\ * ,„x 

Kab = Kab I 1 I + TT a b , (6) 

where p and P are the total matter density and pressure, q a is the total matter momentum density, and Tr a b is the 
matter anisotropic stress. We have neglected terms quadratic in q a and Tr a b since these do not contribute to the 
fluctuations in the effective variables in linear perturbation theory. 

In the background, q* = = ir* b . The modified Fricdmann equations are 

# 2 = ^P° S + U - 4 , (7) 
i 3 a z 

H = -y(p eff + P eff ) + J, (8) 
and the KK energy density behaves like 'dark' radiation: 

p* cx 1 . (9) 
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The (l+3)-covariant description of tensor perturbations in general relativity [13, 14] generalizes naturally to 
brancworld cosmology [15]. Tensor modes on the brane are characterized by the transverse traceless shear a a b, 
which is a tensor potential for the electric and magnetic Weyl tensors E a b and H a b. The sources for linearized tensor 
modes are the matter anisotropic stress ir a b and the KK anisotropic stress ir* b . The Boltzmann equation governs ir a b. 
As pointed out above, the KK anisotropic stress, which carries the imprint of 5D gravitational wave modes on the 
brane, is not determined by 4D equations on the brane. For the moment, we assume that 7r* b is given from a solution 
to the 5D equations. Below we will discuss what this entails and how one can approximate the solution. 

The transverse traceless quantities can be expanded in electric (Q^) an d magnetic (Q^b) parity tensor harmon- 
ics [13], with dimensionless coefficients: 
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H k Q a k b + H k Q^ 
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TTab = P 



<b = PY.{<Qab+**kQa, 
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Using H a b = curl cr a b, we arrive at the coupled equations 
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-pir k = n 2 (2 - 
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^ (E' k + HE k ) - k + ^ - y 7 p) °k + yM - y (37 - WP*k 
^-{2k 1P 2 a k - (37 - 2)pV fc - [3V - (3 7 - 2)(6 7 - l)H]p\ k } 
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where a prime denotes d/dr, with r conformal time, TL = a' /a, and the (non-constant) parameter 7 is defined by 
P = (7 — l)p. Equations (15) and (16), with all braneworld terms on the right-hand sides, determine the tensor 
anisotropies in the CMB, once ir k and 7r^ are given. The former is determined by the Boltzmann equation in the 
usual way [13], except that the background dynamics are altered by braneworld effects. The latter requires a solution 
of the 5D perturbation equations. 

The solution for the non-local anisotropic stress will be of the form 



7^(t)cx J dfg(r,f)F[ 



TTfe, CTfe 



(17) 



where Q is a retarded Green's function evaluated on the brane. The functional F is known in the case of a Minkowski 
background [17], but not in the cosmological case. (An equivalent integro-differential formulation of the problem is 
given in Ref. [18]; see also Ref. [19].) Once Q and F are determined, Eq. (17) can in principle be incorporated into a 
modified version of Boltzmann codes such as CAMB [8] or CMBFAST [20]. It remains a major task of braneworld 
cosmological perturbation theory to find this solution, or its equivalent forms in other formalisms. In the meanwhile, 
in order to make progress towards understanding braneworld signatures on tensor CMB anisotropies, we can consider 
approximations to the solution. 

The nonlocal nature of ir^", as reflected in Eq. (17), is fundamental, but is also the source of the great complexity of 
the problem. The lowest level approximation to 7r£ is local. Despite removing the key aspect of the KK anisotropic 
stress, we can get a feel for its influence on the CMB if we capture at least part of its qualitative properties. The key 
qualitative feature is that inhomogeneity and anisotropy on the brane are a source for KK modes in the bulk which 
"backreact" [15] or "feed back" [16], onto the brane. The transverse traceless part of inhomogeneity and anisotropy 
on the brane is given by the transverse traceless anisotropic stresses in the geometry, i.e. by the matter anisotropic 
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stress Tt a b and the shear anisotropy a a b- The radiation and neutrino anisotropic stresses are in turn sourced by the 
shear to lowest order (neglecting the role of the octupole and higher Lcgcndre moments). 

Thus the simplest local approximation which reflects the essential qualitative feature of the spin-2 KK modes is 



K 2 7T, 



ab 



-{Ha,*, C' = 0, 



(18) 



where £ is a dimensionless KK parameter, with £ — corresponding to no KK effects on the brane, and ( = = A -1 
giving the general relativity limit. [Note that for tensor perturbations, where there is no freedom over the choice of 
frame (i.e. u a ), there is no gauge ambiguity in Eq. (18). However, for scalar or vector perturbations, this relation 
could only hold in one frame, since 7r* b is frame-invariant in linear theory while a a \, is not.] 

The approximation in Eq. (18) has the qualitative form of a shear viscosity, which suggests that KK effects lead to 
a damping of tensor anisotropics. This is indeed consistent with the conversion of part of the zero-mode at Hubble re- 
entry into massive KK modes [5, 7]. The conversion may be understood equivalently as the emission of KK gravitons 
into the bulk, and leads to a loss of energy in the 4D graviton modes on the brane, i.e. to an effective damping. The 
approximation in Eq. (18) therefore also incorporates this key feature qualitatively. 

With this first approximation, we can close the system of equations on the brane by adding the equation 



K 2 pir* k 



k 



(19) 



We will also assume K = = p* in the background. The parameter ( (together with the brane tension A) then 
controls braneworld effects on the tensor CMB anisotropics in this simplest approximation. 



III. 



CMB TENSOR POWER SPECTRA 



Ignoring the photon anisotropic stress (i.e. ir a b = 0), the variable u k = a 1+i */ 2 o k satisfies the equation of motion 

( a -i-C/2)" 



< + 



k z + 2K 



-l-C/2 



Mfe = 0, 



(20) 



where we have used Eq. (18). In flat models (K = 0) on large scales there is a decaying solution a k oc a~( 2+ ^ . 
Since Eq. (20) contains no first derivative term the Wronskian is conserved. On large scales we can use the solution 
a k oc a^^ 2+ ^ to write the conserved Wronskian as W = a' k + (2 + ()TCa k . (The Wronskian vanishes in the decaying 
mode.) Integrating gives the following two independent solutions on large scales in flat models: 



Ok = 



A k a-( 2 +0 , 
B k a-( 2 +0 f dfa(f) 



(21) 



where A k and B k are constants of integration. If we let £ — > 0, we recover the results in Ref. [9] for the general 
relativity case. 

The conserved Wronskian is proportional to the metric perturbation variable, Ht, characterising the amplitude of 
4D gravitational waves. In flat models H T is related to the covariant variables quite generally by 



(22) 



Ignoring photon anisotropic stress, we can eliminate the electric part of the Weyl tensor via the shear propagation 
equation (15) to find kHxk = — cr^, — (( + 2)H(Tk = —W. The fact that Ht is conserved on large scales in flat models 
in the absence of photon anisotropic stress can also be seen directly from its propagation equation, 

H Tk + (2 + QHH Tk + k 2 H Tk = . (23) 

We can solve Eq. (20) on all scales in the high-energy (p A and a oc r 1 / 3 ) and low-energy (p«A and oar) 
radiation-dominated regimes, and during matter-domination (a oc r 2 ). The solutions are 



Ufe(r) 
Hfe(r) 



fcr 
fcr 
fcr 



ciJi (5+c) (fcT) +c 2 ri (5+c) (fcr) 

c 3 Ji (3 +c)( fcT ) + c 4^i( 3 +c)( fcr ) 
c 5 J| +c (fcr) +c 6 y| +( -(fcr) 



(high energy radiation), 
(low energy radiation), 
(matter domination), 



(24) 
(25) 
(26) 
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where Cj are integration constants. The solutions for the electric part of the Weyl tensor can be found from Eq. (15). 
For modes of cosmological interest the wavelength is well outside the Hubble radius at the transition from the high- 
energy regime to the low-energy. It follows that the regular solution (labelled by c\) in the high-energy regime will 
only excite the regular solution (c 3 ) in the low-energy, radiation-dominated era. Performing a series expansion, we 
arrive at the appropriate initial conditions for large-scale modes in the low-energy radiation era: 

H "-'-5^ + sp + O^O +01(trn (27) 



fcr k 3 r 



3 

5l 



- fc = -3Tc + 2( 3 + c)(5 + C) +Q[(fcr)] ' (28) 

(4 + Q (fcr) 2 (8 + Q ,4, (9q] 

Ek ~ 2(3+7) " 4(3 + C)(5 + C) + ° [{kT) 1 (29) 

In the limit ( — > 0, we recover the general relativity results [13]. 

For modes that are super-Hubble at matter-radiation equality (i.e. fcr cq <C 1), the above solution joins smoothly 
onto the regular solution labelled by c 5 in Eq. (26). For fcr cq >• 1, the shear during matter domination takes the form 



= -2^ +c r(§ + C)(fcr)-<i + J5 +C (fcr). (30) 



In the opposite limit, the wavelength is well inside the Hubble radius at matter-radiation equality. The asymptotic 
form of the shear in matter domination is then 

Uk „ F[ ^ (3 + C)] ( ^ (jfer)-(i+C/2) sin(fcT _ nC/4 y (31) 



We use the initial conditions, Eqs. (27)-(29), in a modified version of the CAMB code to obtain the tensor tem- 
perature and polarization power spectra. The temperature and electric polarization spectra are shown in Figs. 1 
and 2 for a scale-invariant initial power spectrum. The normalisation is set by the initial power in the gravity wave 
background. Figs 1 and 2, together with Eqs. (24)~(29), are the main result of this work, and we now discuss the 
physical conclusions following from these results. 



IV. DISCUSSION 



As expected, we find that the power spectra are insensitive to high-energy effects, i.e. effectively independent of 
the brane tension A : the ( = curve in Fig. 1 is indistinguishable from that of the general relativity model (both 
power spectra are identical at the resolution of the plot). For the computations, we have used the lowest value of the 
brane tension A = (100 GeV) 4 , consistent with the tests of Newton's law. 

There are three notable effects visible in Figs. 1 and 2, arising as the physical consequences of our approximate 
model of the KK stress: (i) the power on large scales reduces with increasing KK parameter £; (ii) features in the 
spectrum shift to smaller angular scales with increasing £; and (iii) the power falls off more rapidly on small scales as 
£ increases. Neglecting scattering effects, the shear is the only source of linear tensor anisotropies (see e.g. Ref. [13]). 
For 1 <C I < 60 the dominant modes to contribute to the temperature C;s are those whose wavelengths subtend an 
angle <~ l/l when the shear first peaks (around the time of Hubble crossing). The small suppression in the C;s on large 
scales with increasing £ arises from the reduction in the peak amplitude of the shear at Hubble entry [see Eq. (30)], 
qualitatively interpreted as the loss of energy in the 4D graviton modes to 5D KK modes. 

Increasing £ also has the effect of adding a small positive phase shift to the oscillations in the shear on sub-Hubble 
scales, as shown e.g. by Eq. (31). The delay in the time at which the shear first peaks leads to a small increase 
in the maximum I for which 1(1 + 1)C/ is approximately constant, as is apparent in Fig. 1. The phase shift of the 
subsequent peaks in the shear has the effect of shifting the peaks in the tensor Qs to the right. For I > 60 the main 
contribution to the tensor anisotropies at a given scale is localized near last scattering and comes from modes with 
wavenumber k <~ I /to, where To is the present conformal time. On these scales the gravity waves have already entered 
the Hubble radius at last scattering. Such modes are undergoing adiabatic damping by the expansion and this results 
in the sharp decrease in the anisotropies on small scales. Increasing the KK parameter £ effectively produces more 
adiabatic damping and hence a sharper fall off of power. The transition to a slower fall off in the C;s at I <~ 200 is 
due to the weaker dependence of the amplitude of the shear on wavenumber at last scattering for modes that have 
entered the Hubble radius during radiation domination [21]. [The asymptotic expansion of Eq. (30) gives the shear 



6 



ID' 2 F 1 




2 Q-9 I I I I I 'I' | | | | ill | | | | ill | | | | II 

10 1 10 2 10 3 10 4 

1 



FIG. 1: The temperature power spectrum for tensor perturbations in braneworld models using the approximation in Eq. (18), 
with £ the dimensionless KK parameter. Models are shown with £ = 0.0, 0.1, 0.25, 0.5 and 1.0. The initial tensor power 
spectrum is scale invariant and we have adopted an absolute normalisation to the power in the primordial gravity wave 
background. The background cosmology is the spatially flat ACDM (concordance) model with density parameters tl b — 0.035, 
fl c — 0.315, Q.a = 0.65, no massive neutrinos, and the Hubble constant Ho = 65 km s _1 Mpc _1 . 



amplitude oc k~^ 2+ ^ at fixed r, whereas for modes that were sub-Hubble at matter-radiation equality Eq. (31) gives 
the amplitude oc fc~( 1+ ^/ 2 ).] 

Similar comments apply to the tensor electric polarization Cf, shown in Fig. 2. As with the temperature 
anisotropics, we see the same shifting of features to the right and increase in damping on small scales. Since po- 
larization is only generated at last scattering (except for the feature at very low I that arises from scattering at 
rcionization, with an assumed optical depth tc — 0.03), the large-scale polarization is suppressed, since the shear 
(and hence the temperature quadrupole at last scattering) is small for super-Hubblc modes. In matter domination 
the large-scale shear is <Jk = — fcr/(5 + 2(); the reduction in the magnitude of the shear with increasing KK parameter 
£ is clearly visible in the large-angle polarization. The braneworld modification to the tensor magnetic polarization 
Cp has the same qualitative features as in the electric case. 

In principle, observations can constrain the KK parameter which controls the generation of 5D modes within our 
simplified local approximation, Eq. (18). The other braneworld parameter A, the brane tension, is not constrained 
within our approximation. In practice, the tensor power spectra have not been measured, and the prospect of useful 
data is still some way off. What is more important is the theoretical task of improving on the simplified local 
approximation we have introduced. This approximation has allowed us to encode aspects of the qualitative features 
of braneworld tensor anisotropics, which we expect to survive in modified form within more realistic approximations. 
However, a proper understanding of braneworld effects must incorporate the nonlocal nature of the KK graviton 
modes, as reflected in the general form of Eq. (17). It is also necessary to investigate the scalar anisotropies, which 
have a dominant contribution to the measured power spectra. These may reveal new braneworld imprints that are 
more amenable to observational testing. 
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FIG. 2: The electric polarization power spectrum for tensor perturbations for the same braneworld models as in Fig. 1. 
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